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Topological media are systems whose properties are protected by topology and thus are robust 
to deformations of the system. In topological insulators and superconductors the bulk-surface and 
^ , bulk-vortex correspondence gives rise to the gapless Weyl, Dirac or Majorana fermions on the 

^ ' surface of the system and inside vortex cores. Here we show that in gapless topological media, the 

, bulk-surface and bulk-vortex correspondence is more effective: it produces topologically protected 

gapless fermions without dispersion - the flat band. Fermion zero modes forming the flat band are 
localized on the surface of topological media with protected nodal lines [l|,0| and in the vortex core 
in systems with topologically protected Fermi points (Weyl points) 0| . Flat band has an extremely 
' O ' singular density of states, and we show that this property may give rise in particular to surface 

^ , superconductivity which could exist even at room temperature. 

o 

> ' I. INTRODUCTION 

o 

^\ . Topological matter is characterized by a nontrivial topology in momentum space [^-|7[. This topology may be 
' represented by the momentum-space invariants as depicted in Fig. 1. They are in many respects similar to the 
04 real-space invariants which describe topological defects in condensed matter systems and in particle physics. In 
■ particular, the Fermi surface in metals is topologically stable, because it is analogous to the vortex loop in superfluids 
T— I , or superconductors!^. In the same way, the Fermi point (the Weyl point) corresponds to the real-space point defects, 
^ [ such as hedgehog in ferromagnets or magnetic monopolc in particle physics. The fully gapped topological matter, 
such as topological insulators and fully gapped topological superfluids represent skyrmions in momentum space: they 
have no nodes in their spectrum or any other singularities, and they correspond to non-singular objects in real space 
5^ , - textures or skyrmions (Fig. 1 top right). 

Recently the interest to topological media has been mainly concentrated on the fully gapped topological media, such 
as topological insulators and superfluids or superconductors of the ■^He-B type. These systems contain topologically 
protected gapless fermions on the surface^, Q, and in the core of topological objects some of which have an 
exotic Majorana nature. 

The gapless topological media also exhibit exotic fermion zero modes with interesting properties. In particular 
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they may have a Fermi arc (Fermi surface which terminates on a monopole in Fig. 1 bottom right)) Ill4l3l|. and 



a dispersionless branch of the spectrum with zero energy - the flat band [ll-Q, 14-16[. Historically the flat bands 
were first discussed in relation to Landau levels, but they may emerge even without magnetic fields. They were 
suggested in strongly interacting systems 17-2^, in the core of quantized vortices 21 , in 2+1 dimensional quantum 



22j , in rhombohedral graphite 23 
2M- Flat band 
(Fig. 1 bottom 




field theory which is dual to a gravitational theory in the anti-de Sitter background^ 
and on graphene edge ,24] , and on the surface of superconductors with gap nodes in the bulk [l[ 
is the momentum-space analog of a domain wall (soliton) terminating on a half-quantum vortex . 
left). The topologically protected flat bands, which we discuss here, as well as the Fermi arc add a new twist in the 
investigation of the 3-dimensional topological matter, shifting the interest from the topological insulators and fully 
gapped superfluids/superconductors to their gapless 3-dimensional counterparts, such as superfluid ^He-A, graphite, 
topological semi-metals and gapless topological superconductors. Dispersionless bands may serve as a good starting 
point for obtaining interesting correlated and symmetry breaking states 25 1. 
This paper is based on the earlier arXiv version [2^. 



II. SURFACE FLAT BAND IN A SEMIMETAL 

Consider the semimetal with topologically protected nodes in the form of a spiral as in Figs. 2a and 3a. The 
topological invariant in the bulk, supporting the existence and topological stability of the nodal line and as a result 
of the flat band with respect to interactions, is the contour integral in momentum space 

N2 = -^ tT <f dla,H-'ViH , (2.1) 

where H(jp) is the effective matrix Hamiltonian (inverse Green's function at zero energy), and is pseudospin. The 
nodal line is stable if the integral over the contour C around the line is nonzero. On the other hand one may choose 
the contour C as a straight line along the direction pz normal to the surface, i.e. at fixed momentum p|| = (px,Py) 
along the surface. Due to periodic boundary conditions, the points Pz — '^i^ ja are equivalent and the contour of 
integration C forms a closed loop, giving integer values to the integral iV2(p||) if the integration path does not cross 
the point in the bulk where the energy is zero. This integral 7V2(p||) = 1 for any point p|| within the projection 
Sp of the spiral on the surface, and iV2(p||) = outside this region. The states with momentum py G Sp cannot 
be adiabatically transformed into states in topologically trivial (iV2 = 0) media, and therefore a surface state with 
zero energy is formed for all momenta within . P. B u\ This is in contrast with topological insulators, where this 
type of momentum-space invariants can be defined only for some particular values of momenta (but in principle, the 
insulators with topologically protected surface flat bands are possible, we are indebted to A. Kitaev for this comment). 

Such surface band emerges for example in the system of the rhombohedral stacking (123123. . . ) of graphene 
layers m the limit of a large number M of layers. Simultaneously the nodal line is formed in the bulk, which is 
the source of the topological protection of the surface band. The robustness of the flat band to disorder, represented 
by stacking faults, is demonstrated in Fig. 4c. The flat band remains even if the nodal line is transformed to the Fermi 
surface tube as in rhombohedral graphitej29||, if the invariant as the integral around the Fermi tube is conserved. 

In an analogous way stacking layers of topological insulators leads to the formation of a semimetal with two Weyl 
points in the bulk and a Fermi arc on the surface [l^ • The same happens in the process of the growth of the thickness 
of a ^He-A film, when two Weyl points are formed in the bulk fi^, leading to a Fermi arc on the surface 



III. FLAT BANDS ON THE SURFACE OF CUPRATE SUPERCONDUCTORS AND GRAPHENE 

EDGE 



Dirac points in graphene and nodal lines in high- Tc superconductors have also nontrivial topological charges A'^2 = il 
(see e.g. 1^ ^)- This results in the flat bands on the lateral surface of the superconductor and on the special edges 
of graphene [2J|. Figure 36 demonstrates the ID flat band on the zig-zag edge of graphene and flgure 3c - the 2D flat 
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Fermi surface: vortex ring in p-spa ce fully gapped topological matter: 




Weyl point - hedgehog in p-space 
3He-A, vacuum of SM, topological semimetals 



CO 




flat band (Khodel state): Ji- vortex in p-space 



FIG. 1: Topological matter, represented in terms of topological objects in momentum space, {top left): Fermi surface is the 
momentum-space analog of the vortex line: the phase of the Green's function changes by 27r A'^i around the element of the line in 
(oj, p)-space. (top middle): Fermi point (Weyl point) is the counterpart of a hedgehog and a magnetic monopole. The hedgehog 
in this figure has integer topological charge N3 = +1, and close to this Fermi point the fermionic quasiparticles behave as Weyl 
fermions. Nontrivial topological charges A''i and in terms of Green's functions support the stability of the Fermi surfaces and 
Weyl points with respect to perturbations including interactions 0, [23|. {top right): Topological insulators and fully gapped 
topological superfiuids/superconductors are textures in momentum space: they have no singularities in the Green's function 
and thus no nodes in the energy spectrum in the bulk. This figure shows a skyrmion in the two-dimensional momentum space, 
which characterizes two-dimensional topological insulators exhibiting intrinsic quantum Hall or spin-Hall effect, {bottom left): 
Flat band emerging in strongly interacting systems This dispersionless Fermi band is analogous to a soliton terminated 

by half-quantum vortices: the phase of the Green's function changes by tt around the edge of the fiat band {bottom right): 
Fermi arc on the surface of ^He-A and of topological semi-metals with Weyl points serves as the momentum-space 

analog of a Dirac string terminating on a monopole. The Fermi surface formed by the surface bound states terminates on the 
points where the spectrum of zero energy states merge with the continuous spectrum in the bulk, i.e. with the Weyl points. 



band on commensurate surface (1,1,0) of d-wave superconductor. Projections of the nodes on the surface determine the 
boundaries of the regions at which the integral A'^2(P||) takes a definite value: 0, -1-1 or —1 (here p|| = Py is momentum 
along the graphene edge and py = {px,Pz) for momentum along the surface of the superconductor; integration in (|2.ip 
is now along the momentum p±^ perpendicular to the surface/edge). For any point p|| for which Af2(P||) 7^ 0, there is 
a surface bound state with exactly zero energy e(p||) = 0, i.e. flat band. For cuprate superconductors, this flat band 
corresponds to the infinite set of Andreev bound states obtained in the semiclassical approximation and observed in 
grain-boundary junctions, see e.g. [32h34| . The topological analysis demonstrates that even beyond the semiclassical 
approximation the Andreev bound states form the flat band with exactly zero energy in some regions of momentum 
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FIG. 2: Topological media which exhibit topologically protected dispersionless fermionic spectrum with exactly zero energy - 
the flat band, (a): Three-dimensional topological semi- metal. Flat band can be obtained by stacking of graphene layers, each 
alone represented by Dirac points. In case of proper stacking, represented by the coupling of the pseudospins in each layer, the 
nodal line is formed in the bulk in the form of a spiral which ends up on the faces of the Brillouin zoneQjEsi- This nodal line has 
a non-zero topological charge ~ 1. This nontrivial charge protects the surface states with zero energy e{px,Py) — 0, in the 
whole region within the projection of the spiral on the surface, see Fig. 3a. (6): Dirac points in graphene {top: real-space crystal 
lattice; bottom: reciprocal space). They lead to the flat band on the zig-zag edge, see Fig. 3b. (c): Cuprate superconductors 
have also topologically protected nodal lines [l^. [iH. They lead to the fiat band on a side surface, see Fig. 3c. (d): Superfluid 
^He-A has two three-dimensional Weyl points, with Ns = +1 and N3 = —1. The projections of the Weyl points on the direction 
of the vortex line determine the boundaries of the region where the spectrum of fermions bound to the vortex core is exactly 
zero, e(pz) = 0, see Fig. 3d and Fig. 5. 



flat band: half-quantum vortex in p-spare flat band m the vortex core 




FIG. 3: Flat band formation, (a, top): Two-dimensional flat band appears on the surface of gapless systems with a topologically 
protected nodal line 1, 2]. (a, bottom): The nodal spiral in a semi-metal. This nodal line has a non-zero topological charge in 
momentum space, ~ 1, and this charge protects the surface states with zero energy e{px,Py) = 0, in the whole region within 
the projection of the spiral on the surface. (6): One-dimensional topologically protected flat band emerges on the zig-zag edge 
of graphene. Projections of the Dirac points on the edge determine the boundaries of the flat band (c): Nodes (nodal lines) 
in cuprate superconductors give rise to the two-dimensional flat bands on the lateral surface [i^. Projections of the nodes on 
the surface determine the boundaries of the flat bands, (d): One-dimensional topologically protected flat band emerges in the 
core of linear topological defects (such as the vortex with winding number n = 1 in real space in the bottom right comer) in 
three-dimensional topological matter with Weyl points. The projections of the Weyl points on the direction of the vortex line 
(along the z-axis) determine the boundaries of the region where the spectrum of fermions bound to the vortex core is exactly 
zero, e(pz) = 0, see also Fig.5. 
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FIG. 4: Formation of the flat band and nodal line by stacking M graphene-like layers in the limit of large M . Flat band is 
formed from the multiple Dirac point on the surface, while the nodal line is formed from the gapped states in the bulk, (a): 
Formation of zero energy states in the bulk. Minima of the bulk energy spectrum at fixed P|| = (px,Py) approach zero when 
M — >■ CO. These zeroes correspond to the projection of the spiral in Fig. 3a onto the surface, (b): In the same limit Af — )■ oo the 
surface flat band is formed. Its boundary coincides with the projection of the nodal spiral onto the surface. This demonstrates 
the bulk-surface correspondence for the flat band, (c): Robustness of the surface flat band to the stacking fault disorder. The 
blue curves show the flat band eigenvalues vs. px for py = in the case when there is only one coupling t+ = t in the model 
The red curves show the eigenvalues in the case of Gaussian-distributed disorder in both t+ and couplings with a standard 
deviation equal to 0.5i and the black curves in the case when the disorder in both couplings equals the average coupling t. In 
the previous case the flat band survives, whereas in the latter case of strong disorder it is destroyed. The results are obtained 
for Af = 200 layers and averaged over disorder. 



space if the surface has an orientation commensurate with the crystal. In this case the flat band consists of Majorana 
fermions, and for the considered (1,1,0) surface there are boundaries which separate flat bands with different topology, 
N2 = +1 and N2 = -1. 



IV. FLAT BAND IN A VORTEX CORE 



The topologically protected dispersionless spectrum with zero energy appears also in linear topological defects - 
vortices. The topological protection is provided by the nontrivial topology of three-dimensional Weyl points in the 
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FIG. 5: Schematic illustration of the spectrum of bound states E{pz) in the vortex core. The branches of bound states terminate 
at points where their spectrum merges with the continuous spectrum in the bulk. The flat band terminates at points where 
the spectrum has zeroes in the bulk, i.e. when it merges with Weyl points. It is another p-space analog of a Dirac string 
terminating on a monopole, see Fig. 1 bottom right for the Fermi arc. 

bulk. This bulk-vortex correspondence is illustrated in Fig. 3d for the case of a pair of Weyl points with opposite 
topological charges = ±1. The projections of these Weyl points on the direction of the vortex line determine the 
boundaries of the region where the spectrum of fermions bound to the vortex core is exactly zero, E{p^) = 0. Such 
flat band first obtained in Ref. [21] for the noninteracting model is not destroyed by interactions. The spectrum of 
bound states in a singly quantized vortex in ^He-A is illustrated in Fig. [5] 

V. DENSITY OF STATES, ENTROPY AND SPECIFIC HEAT OF A FLAT BAND 

The important consequence of the flat band is the anomalously singular density of states (DOS) of electrons. The 
formation of a singular DOS with growing the number of layers M can be found in the model 2] in which the spectrum 
of electrons in the surface layer contains the multiple Dirac point with energy e(p||) = i(|p|| |/po)*^: where t is the 
intcrlayer coupling. For M > 2 this system has a singular DOS z^(e) ^ g^/Af-i ^^j. j^j ^ 2 it is ^{e) = l/(27rAfe) 
at e = 0. Assuming for simplicity the spectrum e(P||) ~ ^(|P|| |/Po)^^ for all P < Po one obtains the following entropy 
of the surface flat band for large M 

S^-^%(l^—lni]ln2 , l»Af»21n-, (5.1) 

where Sp = irp^ is the area of the flat band in momentum space, and A the area of the surface in real space. The first 
term describing the flat band limit can be understood as the entropy related to occupying or not occupying the phase 
space volume ASp/ {2'Kfif' of zero energy states. The second term leads to a finite heat capacity at low temperatures, 
Cm{T) — >■ ^45^ ln2/(27r^fi,^M), which in the range of (|5.ip dominates over the bulk contribution. The exact spectrum 
gives the same result which corresponds to the equipartition law with the energy nn2/(2M) per each degree of 
freedom in the flat band. 

VI. SUPERCONDUCTIVITY 

Singular DOS gives rise to superconductivity with an enhanced transition temperature. Already for M = 4, 
i'(e) ^ corresponds to the DOS at the edge of a one-dimensional band, which for a long time has been discussed 
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as a possible source of enhancement of the superconducting transition temperature, see e.g. [35|. For general M the 
self-consistency equation for the gap on the surface reads 



For simplicity we chose the pairing potential Vp = g9{pm — p) with pm ^ po- For the flat band with e^(P||) = one 
obtains 

A = ^tanh— . (6.2) 

For T = we get A(0) = gp^/SiTh^ and the corresponding transition temperature is Tc — A(0)/2. For finite M one 
obtains a power-law behavior of the superconducting gap at T = 0, 

M A/ -1-2 



For M ^ 1 it approaches the linear law in Eq. (|6.2p . For cut-off pm > Po the more detailed calculation gives 
A(0) = gSp/lGir^h? , where again Sp is the area of flat band in momentum space. The linear dependence of the 
gap and of the critical temperature on the coupling should be contrasted with the exponential suppression of Tc in 



conventional bulk superconductivity, Tc ~ exp(— l/z^(0)(7). The details can be found in 15l. Il6|. 

The flat band should also emerge at the twin boundaries, if the neighboring twins have opposite topological polar- 
ities, i.e. opposite helicities of the nodal spiral. Such flat bands also support superconductivity, and thus a proper 
arrangement of many twins in the bulk would give bulk superconductivity with a large transition temperature. Ex- 
istence of localized superconducting domains at elevated temperatures has been suggested for pure graphite and for 
graphite-sulfur composites 36, S^], which may be related to the topology of nodes discussed here. 

The strong dependence between the superconducting critical temperature and the coupling energy due to the 
singular density of states in the flat bands may open the route for room temperature superconductivity. However, 
also other types of broken symmetry states 38|, such as magnetic instabilities 39 1 or exciton BEc[4^ are possible. 



depending on the nature and magnitude of the corresponding coupling constants. 

VII. CONCLUSION 

Topological matter has generic behavior, depending on the topological class. Fermi surface, being topologically 
protected from interaction, results in the robustness of all the generic properties of metals at low temperatures. 
Systems with Fermi points (Weyl points in 3D and Dirac points in 2D) represent other important classes of topological 



matter. Systems with Weyl points (superfluid '^He-A ^ and topological semimetals 12, 13, 41, 42 1) experience 
emergent chiral or Majorana fermions and in addition the effective gauge fields and gravity as collective or artificial 
fields related to the deformations of the Weyl point. There are many reasons to believe that the quantum vacuum 
belongs to the class of topological media with Fermi points. Quantum Hall effect, superfluid ^He-B and topological 
insulators have provided a new class of topological matter, which is characterized by topologically protected surface 
states and quantization of physical parameters, such as Hall and spin-Hall conductivity. The highly degenerate 
topologically protected state - the flat band - is also a generic phenomenon. This new class of topological phenomena 
is waiting for its exploration. 
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